Accurate description of bulk and interfacial properties in colloid-polymer mixtures 
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Large-scale Monte Carlo simulations of a phase-separating colloid-polymer mixture are performed 
and compared to recent experiments. The approach is based on effective interaction potentials 
in which the central monomers of self-avoiding polymer chains are used as effective coordinates. 
By incorporating polymer nonideality together with soft colloid-polymer repulsion, the predicted 
binodal is in excellent agreement with recent experiments. In addition, the interfacial tension as 
well as the capillary length are in quantitative agreement with experimental results obtained at a 
number of points in the phase-coexistence region, without the use of any fit parameters. 

PACS numbers: 82.70.Dd, 61.20.Ja, 82.70.-y 



In addition to commercial applications, mixtures of 
colloids and nonadsorbing polymer are interesting be- 
cause of their analogy to atomic systems p|. Much ef- 
fort has been devoted to understand the phase behavior 
of such mixtures. Of particular interest is phase sepa- 
ration, which occurs when the polymer density and di- 
ameter of gyration are sufficiently large, leading to the 
formation of two coexisting phases: one phase lean in col- 
loids and dense in polymers (the colloidal vapor) and one 
phase dense in colloids and lean in polymers (the colloidal 
liquid). As was shown by Asakura and Oosawa (AO), 
phase separation in colloid-polymer mixtures (CPM) is 
driven by entropy 0]. In the AO description, colloids 
and polymers are treated as effective spheres, assum- 
ing hard-sphere interactions between colloid-colloid and 
colloid-polymer pairs, while the polymers can interpene- 
trate freely. A major advance has been the development 
of a geometry-based density functional for the AO-model 
, which has led to a host of novel and intriguing predic- 
tions regarding interfacial phenomena within this model 
0. However, when comparing to actual experiments, 
quantitative discrepancies arise. An important deviation 
of the AO-model is that it underestimates the polymer 
concentration in the colloidal liquid 5, 6J. This effect is 
more pronounced within the free-volume approximation 
, which is also the bulk limit of the density functional 
of Ref. 0, and it persists when the AO-binodals are ob- 
tained by means of computer simulations 0, Q . 

While the colloid-colloid interaction in realistic sys- 
tems is indeed well described by the hard sphere po- 
tential 9], the colloid-polymer and polymer-polymer in- 
teractions are more complex. To circumvent the short- 
comings of the AO-model, numerous different approaches 
have been employed, both at the effective |l0L UlL Il3| 
and at the monomer-resolved 01111 levels. As a general 
trend, inclusion of polymer nonideality does improve on 
the major drawback of the AO-model, i.e., it yields higher 
polymer concentrations in the colloidal liquid fill 1 1 .' j|] . A 
remarkably accurate method to capture polymer nonide- 



ality is the "polymer as soft colloid" approach [li|. Here, 
the polymers' centers of mass are chosen as effective co- 
ordinates while all fluctuating monomers are canonically 
integrated out. The sought-for effective potentials are ob- 
tained by inverting the correlation functions obtained in 
Monte Carlo (MC) simulations of self-avoiding random 
walks (SAW) on a lattice [H IH El- The pioneering 
work of Bolhuis et al. 01 nas led to the most accurate 
determination of the phase behavior of CPM to date, 
demonstrated by direct comparison to experiments p| . 

Despite the progress in calculating bulk phase dia- 
grams, accurate predictions of the interfacial tension 
7 between coexisting phases, a quantity that plays a 
key role in wetting and interfacial phenomena, remain 
elusive. Theoretical approaches include the square- 
gradient approximation 0, 0, H3, density- functional 
theory [H El H2 and simulations H |H Exper- 
imental data on 7 are hard to obtain, mainly due to 
the very small value of this quantity (7 ~ 1/iN/m). 
In several cases 0, 0, 0| , theoretical predictions have 
been compared to the experimental results of Refs. p4| 
and |25| . These comparisons are carried out by plot- 
ting the theoretical interfacial tension as a function of 
the colloid density gap across the binodal, a procedure 
that tends to obscure the fact that the theoretical and 
experimental binodals can be in considerable disagree- 
ment 0. Recently, Moncho-Jorda et al. [22j employed 
density-functional theory to calculate the interfacial ten- 
sion employing the interactions of Ref. [5j- However, 
they adopted a depletion picture, losing thereby the ef- 
fect of polymer-induced many-body interactions between 
the colloids. Thus, a full two-component treatment of in- 
teracting CPM is necessary in order to capture bulk and 
interfacial behavior quantitatively . 

In this Letter, we demonstrate that simulations of 
CPM using accurate effective interactions predict bulk 
and interfacial properties correctly. We consider the so- 
called "colloid limit" , where the polymer diameter of gy- 
ration <tg is smaller than the colloid diameter a c . Re- 
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cently, experimental measurements of the interfacial ten- 
sion and capillary length in the colloid limit became avail- 
able, to which we can compare 26]. At the same time, 
accurate effective interactions for the colloid limit exist, 
obtained in off-lattice molecular dynamics (MD) simula- 
tions of self-avoiding polymer chains (22) . We choose here 
the central monomer to represent the chain, canonically 
tracing out the remaining, fluctuating monomers. If the 
coarse-graining procedure is carried out accurately, the 
bulk thermodynamics of the mixture should be strictly 
independent of the choice of the effective coordinates [28| . 
Moreover, for athermal solvents and in the limit of long 
chains, details of the microscopic monomer-monomer in- 



teractions become irrelevant |29|. Thus, the phase be- 
havior should be independent of whether one adopts for 
the polymers a lattice, SAW- model as done in Ref. p| , or 
the approach at hand. We will explicitly check whether 
this requirement is fulfilled in what follows. 

The colloids are represented by their centers and r 
denotes the distance between any effective coordinates 
in the corresponding effective interactions Vij{r), with 
i, j = c,p. The colloid-colloid interaction V cc (r) is given 
by a hard-sphere potential of diameter a c . The colloid- 
polymer interaction, V cp (r) diverges for r < cr c /2 and for 
larger separations it reads as |27j : 
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where a p a typical length scale given by a p = 0.66ctg, 
(3 = (fceT) -1 , T the temperature, fee the Boltzmann 
constant, £i = 1/(1 + 2K 2 a 2 ), = (V^£i/ Kcr p)[ 1 - 
erf(/t(7 p )] exp(K 2 gp , and parameters A and k, deter- 
mined in Ref. 27] by fitting to simulation results. The 
corresponding polymer-polymer interaction is given by: 
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with r obtained by requiring that the effective interaction 
correctly reproduces the second virial coefficient of dilute 
poly mer solutions and resulting in the value t<j p = 1.03 
|27|. To match the experiment of Ref. j2(J , we consider 
a mixture of colloids and polymers with size ratio q = 
cg/cc = 0.56. The parameters are A = 0.46, and Ka p = 
0.52715 H3. 

The binodal and the interfacial tension are obtained 
in the grand canonical ensemble, by MC simulation of 
a mixture of N c colloids and Ap polymers, interacting 
via the above pair potentials. In this ensemble, the 
temperature, the volume V, and the respective fugac- 
ities, z c and z p , of colloids and polymers, are fixed, 
while the number of particles in the system fluctuates. 
We also introduce the colloid and polymer packing frac- 
tions r] CtP = (ira^ G /6)N c , p /V . Since the interactions are 
athermal, the temperature plays no role and the phase 
behavior is set by q and the fugacities. The polymer fu- 
gacity is used as control parameter, analogous to inverse 
temperature in fluid- vapor transitions in atomic systems. 
For a given z p , we measure the distribution P(r] c ), de- 
fined as the probability of observing a system with colloid 
packing fraction rj c . For z p sufficiently far away from the 
critical point, phase coexistence is obtained by tuning z c 
so that P{rj c ) becomes bimodal, with two peaks of equal 
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FIG. 1: Colloid-polymer binodals for q = 0.56 obtained in 
simulations using two different models (ESP and AO). Open 
circles are experimental state-points at which phase separa- 
tion was observed; black triangles are experimental state- 
points at which only one phase was observed. The experimen- 
tal data were taken from Refs. [2(J and |33|. The free- volume 
AO-binodal is also shown. 



area. The peak at low n c corresponds to the colloidal va- 
por phase, the peak at high rj c to the colloidal liquid, and 
the region in between to phase-separated states [sj. The 
average peak locations yield the colloid packing fractions 
in the two phases. The interfacial tension is obtained 
from the average height of the peaks 0, 0HJ. In order 
to simulate efficiently, a grand canonical cluster move is 
used , in combination with a reweighting scheme j^] . 
To obtain the binodal, we vary z p and record the cor- 
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FIG. 2: Interfacial tension 7 against the polymer fugacity 
from our ESP-simulations. The line is a guide to the eye. 



responding densities of colloids and polymers in the two 
coexisting phases. This yields the phase diagram in sys- 
tem representation, which may directly be compared to 
experiments, see Fig. ^ The solid curve is the binodal 
of the present work incorporating the effective soft po- 
tentials (ESP) V cp {r) and V pp (r) above. For compari- 
son, the dashed line shows the binodal of the AO model 
with q = 0.56, obtained following the same grand canon- 
ical simulation procedure. The free volume result is also 
shown. As seen in Fig. ^ ESP interactions give an ac- 
curate description of the experimental binodal. Note in 
particular the significant increase in the polymer density 
at high colloid density in comparison to the AO-result. 

Next, we consider the interfacial tension, and compare 
to the recent experiment of Ref. (2|| . The colloids are 
PMMA spheres (o~ c — 50 nm), mixed with polymer with 
diameter of gyration ctg = 28 nm, and dissolved in de- 
calin at T = 298 K. In Fig.0 we show simulation results 
of the interfacial tension as function of z p . As expected, 
the tension decreases markedly upon lowering z p , van- 
ishing at the critical point (the critical polymer fugacity 
is approximately z PjCr w 3.2). To enable the comparison 
to experiment, we show in Fig. [21 several binodal ticlincs 
obtained in the simulation. At point X, which is close 
to the tieline corresponding to z p = 4.25, the experimen- 
tal interfacial tension equals 7 = 0.16 — 0.2 /xN/m [26j . 
The corresponding interfacial tension in the simulation 
reads 7 = 0.22 /xN/m, which exceeds the experiment by 
only 10%. In contrast, the interfacial tension obtained 
at the same state-point in a simulation of the AO model 
is 7 = 0.07 /iN/m, which underestimates the experiment 
by over 50%. This is due to subtle differences in the 
location and range of the critical regions of the two mod- 
els. More precisely, defining the distance from the criti- 
cal point as t = z p /z PtCt — 1, we obtain for point X using 
ESP interactions t w 0.33, but only t « 0.05 using the 



FIG. 3: Close-up of the vapor branch of the binodal. Open 
squares with tielines (dashed) are binodal points obtained us- 
ing ESP interactions. The polymer fugacity corresponding to 
the tielines reads z p = 7.5 — > 7 — > 6.5 — > 6 — > 5.5 — > 5 — > 
4.5 -> 4.25 -> 4 -> 3.75 -> 3.5 -> 3.45 (from top to bottom). 
The solid curve is the binodal of the AO model. The X-symbol 
(r) c = 0.076, ri P = 0.50) marks a state-point at which the in- 
terfacial tension was measured experimentally. The closed 
circles on the dilution lines A and B are experimental state- 
points at which the capillary length was measured. 



AO model. For the AO model, point X is thus much 
closer to criticality and hence the interfacial tension is 
lower. At the same time, if one calculates the interfacial 
tension for the AO-model using the free volume tieline 
with the same colloid density gap as in the experiment, 
the value 7 = 0.5/zN/m is obtained [2(|. These discrep- 
ancies within the AO-model demonstrate the large ef- 
fect that inaccuracies in the binodals have on 7, as well 
as ambiguities that arise by comparing interfacial ten- 
sions at "rescaled" state-points. In our approach, on the 
contrary, we provide a comparison to the experimentally 
measured interfacial tension in absolute terms, i.e., at the 
same state-point (r] c ,r] p ). 

Finally, we compare to experimental measurements of 
the capillary length l c — ^^/(gAg), with g — 9.81 m/s 2 
the gravitational acceleration, interfacial tension 7, and 
Ag the mass density difference between the colloidal liq- 
uid and vapor phase. For two coexisting phases with 
colloid- and polymer-packing fraction gaps Ar/ C and A?7 p 
respectively, it holds Ag — An c (gc — Qd) + Ar) p m p /!)p, 
with the mass densities g c = 1170 kg/m 3 for PMMA 
and g d — 890 kg/m 3 for decalin, m p = 3.87 x 10~ 22 kg 
the single polymer mass, and v p — ttoq/6 the effective 
single polymer volume |2(j . In Fig. we plot the capil- 
lary length as obtained in the simulation as function of 
z p . The experimental data are also shown, where the 
conversion to z p was performed with the aid of Fig. |3| 
(experimental measurements at z p > 7.5 were converted 
using linear extrapolation). As 7 = gAgl^ and since 
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FIG. 4: Capillary length as function of the polymer fugacity 
obtained in simulations (ESP), as well as in experiments along 
dilution lines A and B of Fig. The line connecting the 
simulation data serves as a guide to the eye. 



Ag is given accurately by our approach (as shown by 
the good agreement with the experimental binodal), the 
agreement with experiment regarding l c directly implies 
agreement with the interfacial tension 7 at all considered 
state-points. Another striking feature of Fig. 0] is the 
remarkably good agreement between simulation and ex- 



periment close to the critical point. For 3D Ising systems, 
the capillary length is expected to vanish at criticality as 
l c ~ t x with the critical exponent A = v — (3/2 w 0.47. 
While the present simulation and experimental data are 
not accurate enough to extract the exponent, the curva- 
ture of the data is certainly compatible with the antici- 
pated exponent. In particular, the data seem to approach 
the z p axis with perpendicular slope. 

In summary, we have demonstrated that the "polymer 
as soft colloid" approach 0], using accurate effective 
colloid-polymer and polymer-polymer interactions, not 
only reproduces the experimental binodal, but also the 
interfacial tension and the capillary length. Note that ex- 
cellent agreement with the experimental binodal was also 
obtained in Ref . |5| , in which the polymer centers-of-mass 
were employed as effective coordinates. Both the effec- 
tive description of Ref. and the one adopted here, 
thus reproduce the correct thermodynamics, providing 
a strong confirmation of the power and self-consistency 
of coarse-graining techniques. We anticipate that a full, 
two-component calculation of the interfacial tension us- 
ing the effective interactions of Ref. [5j will also capture 
the interfacial properties correctly; the latter should be 
the subject of further investigations. Additional experi- 
mental work in measuring interfacial tensions in colloid- 
polymer mixtures is also highly desirable. 
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